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Abstract. We investigate the dependence of spin squeezing on the polar angle of 
the initial coherent spin state j^oi^o) in ^ generalized one-axis twisting model, where 
the detuning 6 is taken into account. We show explicitly that regardless of S and ^o, 
previous results of the ideal one-axis twisting is recovered as long as 6q = it/ 2. For 
a small departure of 6o from 7r/2, however, the achievable variance (y-)inin ^ iV^/'^, 
larger than the ideal case TV^/'^. We also find that the maximal-squeezing time tmin 
scales as N~^/^. Analytic expressions of {V-)rain and tmin arc presented, which agree 
with numerical simulations. 
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1. Introduction 

Spin squeezing arising from quantum correlation of collective spin systems [1], has 
potential applications in high-precision measurement [2], [3] and quantum information 
processes [11 [5], El [3, [HI [9] . Kitagawa and Ueda have studied the squeezing generated by 
a nonlinear Hamiltonian due to the one-axis twisting (OAT) [1]. Starting from a 
coherent spin state (CSS) [lO]: \9o = 7r/2, 0o = 0), the system evolves into spin squeezed 
state (SSS), which shows the reduced variance {V-) below standard quantum limit 
(SQL) — N/A, where is total particle number. The smallest variance (V^_)min ~ N^^^ 
is obtainable at the time scaled as x^min ~ N'"^/'^. 

Possible realization of the OAT-induced squeezing in a two-mode Bose-Einstein 
Condensates (BECs) has been proposed [1], where the self-interaction parameter 
X ~ (oaa + cibb — 2aafe)/2 is inherently aroused from atomic intra- and inter-species 
collisions. Atomic collisions lead to both the squeezing and phase diffusion pTl [T2] . 
The dephasing process destroys phase coherence of the two-component BECs, and thus 
sets a limit to the applications of the condensates in high-precision measurement and 
quantum information processing. A straightforward way to suppress the diffusion is 
the preparation of number-squeezed state, a special case of the SSS with the reduced 
variance along the component. Such a kind of squeezed states have been investigated 
both experimentally [131 [H [IHl [161 [HI [TBI [19] and theoretically [20l [2ll [22l [231 [21] . 

Besides the above schemes that rely on nonlinear interactions of the ultracold atoms, 
spin squeezing can be generated via light-matter interactions [21 [3l [251 [SSI [27] and 
quantum nondemolition measurement [281 [221 [301 [SB [321 [331 IM]- Recently, the OAT- 
induced squeezing has been demonstrated in an ensemble of cesium atoms [311 [32] and 
ytterbium atoms [331 [M]. In their experiments, the CSS with 9q = it/2 was adopted as 
the input state, which is the optimal initial state to obtain the strongest squeezing. Via 
optical pumping, it was shown that 98% atoms are in the CSS [32] . 

In this paper, we investigate the degree of the OAT-induced squeezing for 9q slightly 
departure from 7r/2. A generalized one-axis twisting model: H = 5Jz+xJl considered, 
which is the most important prototype in studying spin squeezing [H H] and quantum 
metrology [351 [26] . We prove exphcitly that without particle losses, the detuning 5 and 
the azimuth angle 0o give vanishing contribution to the squeezing parameter, and the 
ideal OAT-induced spin squeezing can be reproduced as long as = ^/2- As the main 
result of our paper, we investigate the dependence of the variance (l^_)min and the time 
tjnin on the particle number and the polar angle ^o- Our results show that even for 
a small departure of 9q from 7r/2, power rule of the smallest variance (\^_)min changes 
from N^/^ to N'^/^ with the increase of particle number A^. The maximal squeezing is 
achievable at the time that scaled as x^min ~ N~^^^. 

Our paper is organized as follows. In Sec. II, we present general formulas of spin 
squeezing for arbitrary spin-1/2 system. In Sec. Ill, we study quantum dynamics 
of the OAT model, which is exactly solvable for any initial CSS. In Sec. IV, we 
present short-time solutions of the first- and second-order moments of the spin operators. 
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Approximated expression of the reduced variance V- is presented to obtain power rules 
of the maximal squeezing and its time scale tmin- Finally, a summary of our paper is 
presented. 

2. Some formulas of the spin squeezing 

Assume that an ensemble of two-level atoms (i.e., spin 1/2 particles) with ground state 
I a) and excited state \b) can be described by collective spin operator J = Ylk=i 1'^'''^^ 
where a^^^ is the Pauli operator of the kth atom. Spin components of J obey SU(2) 
algebra, [JmyJn^] = iJns for any three orthogonal vectors ni, n2, na. The associated 
uncertainty relation reads (AJn^)^(A Jn^)^ > \\{Jn^)\'^, where the variance is defined as 
usual, (AA)^ = — (\E'|y4|\E')^ for any spin state and operator A. Considering 

the mean spin (J) = {{Jx)i {Jy)i {Jz)), "we choose the orthogonal vectors as 

ni = (— sin 0, cos</), 0) , 

112= (— COS ^ COS 0, — cos6' sin0, sin^) , (1) 
n3 = (sin 6* cos 0, sin 6' sin 0, cos 6*) , 

where the azimuth angles (p = ta.n~^[{Jy) / (J^)], and the polar angle 9 = tan~^[r/( J^)] 
with r = |(J+)| = {{Jx)"^ + {Jy)'^Y^'^ [see Fig. [U^a)]. For arbitrary spin state it is 
easy to prove that the mean spin (J) is along the direction, with the length of the 
mean spin R= |(J)| = (Jns) [see Append.]. Now, let us consider the CSS [TO] : 

1^,0) = e-'''-^\j,j) = e^e(^^-^-^.-«^)|j,j), (2) 

which is eigenstate of Jns with eigenvalue j = N/2 (where A^ is total particle number), 
and thus {Jns) = I (J) I = j- In single-particle picture, the CSS can be rewritten as a 
direct product, 16*, 0) = nfcLi[cos(^/2)|6)fc -|- e*'''' sin(6'/2)|a)fc], where \a)k and \b)k are 
ground and excited states of the kth atom. Such a quantum uncorrelated state obeys 
the minimal uncertainty relationship: (AJ^)^ = (AJ^^)^ = ^K-^ns)! = j/2, where the 
value j/2 is termed as the SQL. 

Since the mean spin is parallel with n^, one can introduce any spin component 
normal to the mean spin as 

= J ■ n^ = Jn^ COS ip + sin ip, (3) 

where the unit vector = nicos^/' -|- n2sin'?/', with ip, being arbitrary angle between 
ni and n^,. For any spin state |\&), we have (J^) = and therefore, the variance of 
reads 

( A J^)2 = ^[C + A cos{2tlj) + B sm{24j)] , (4) 

where the coefficients A = {J^^ — J^^), B = {Jn^Jn-z + Jn^Jni), and C = (J^^ + J^^) = 
j{j + 1) — (Jns)- Another orthogonal spin component with respect to and its variance 
can be obtained by replacing ip with ip + 7i/2. For the CSS |6',0), it is easy to verify 
that the coefficients A = B = and C = j, which gives the variance (AJ^)^ = j/2, 
indicating isotropically distributed variances of the CSS P, as shown in Fig. [l](b). 
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(a) (b) (c) 

Figure 1. (color online) Husimi Q function: Q{d,(f>;t) = \{9,4>\'^{t))\'^ on the Bloch 
sphere for j — 30. (a) the unit vectors ni (green). n2 (blue), (red), as defined in 
Eq. (1). (b) the initial CSS 16*0 = tt/S, 0o = 7r/3). (c) the SSS generated by the OAT 
Haniiltonian H — xJz time imin = 0.043x~^, where x^min is the time scale to attain 
the strongest squeezing. For large j, it is given by Eq. (j28p . 



A spin-squeezed state (SSS) is defined if tfie variance of one spin component normal 
to the mean spin is smaller than the SQL [1], i.e., (AJ^)^ < j/2. The SSS has anisotropic 
variances distribution in a plane normal the mean spin [see Fig. [11(c)]. Optimally 
squeezed angle ipop is obtained via minimizing (AJ^)^ with respect to ip, yielding 
tan(2V'op) = B/A, so cos(2V'op) = ±A/VA'^ + B"^ and sin(2?/'op) 



Substituting these results into Eq. 
[3 [221 [23] 

1 

2 



C±VA^TW 



we obtain the reduced and the increased variances 



(5) 



where the reduced variance = (AJ^)^ corresponds to the squeezing along with 
= 'ipop = [tt + tan~"^(i3/^)]/2; while the increased variance 1+ gives the so-called anti- 
squeezing for the angle ip = i^op + 7r/2. The degree of spin squeezing can be quantified 
by the normalized variance 

C 



2V^ 



(6) 



J J 

For the CSS, the variances = V+ = 3/2 and = i; while for the SSS, < 1. it 
should be mentioned that the coefficients A, B, and C depend only on five quantities 
[see Append. A]: (J^), (J+), {J^), {J+)^ and (J+(2J^ + 1)), from which one can solve the 
mean spin (J) and the squeezing parameter In addition, there are several definitions 
of the squeezing parameter. According to Wineland et al. [2], the squeezing parameter 
is defined as 



2j 



3' 



(7) 



l(J)P ' l(J)P 

which closely relates to both frequency resolution in spectroscopy [2] and many-body 
quantum entanglement [1]. 
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3. Generalized one-cLxis twisting model and its exact solutions 

The above formulas are valid for any spin-1/2 system with SU(2) symmetry. As an 
example, we consider a two-component BECs [371 EH] confined in a deep 3D harmonic 
potential. The total system can be described by the two- mode Hamiltonian {h = 1) 
[39]: 



H = + u;,N, + U,,NaN, + ^{a)f{af + (8) 

where a, 6, and Ni {i = a, b) are the annihilation and number operators for the 
two internal states \a) and \b), Ui are single-particle kinetic energies, and Uij = 
(Airaij/M) J d^r\^Q{r)\'^ are atom-atom interaction strengthes. For a conserved total 
particle number N = Na + Nf,, the two-mode model can be rewritten as if = 6Jz + xJ^i 
where the detuning 5 = u^- u}a + {Ubb - Uaa){N - l)/2, and x = {Uaa + Ubb - 2Uab)/'2. 
Angular momentum operators J+ = {J-Y = u^a and = {Nb — Na)/2, satisfying SU(2) 
algebra. 

Assumed that the two-mode system evolves from the CSS, 1^1/(0)) = 16*0,00) = 
Cm(0)|j, m), with the probability amplitudes [10] 



{j+m)\{j-m)\ \2j \2) 

where the polar angles ^o and 0o determine population imbalance and the relative phase 
between the two internal states [iQl HI] . The state vector at any time t reads 

l*W)= E c™e-^^'"^+>^'"'^*|j,m), (10) 

m=-j 

where the self-interaction x scrambles phase of each number state \j, m) , and leads 
to spin squeezing [Tj H] and phase diffusion [TT] of the two-mode BEC. In theory, 
the diffusion is quantified by correlation function {o^d) (i.e., which decays 

exponentially with the time scale x^d = for 9q = 7i/2. Such a kind of the dephasing 
process has been observed in experiment by extracting the visibility of the Ramsey fringe 

As an ideal case, spin squeezing induced by the OAT Hamiltonian xJz has been 
investigated for the initial CSS \6o = 7r/2,0o = 0) [1]. For this special CSS, it was 
shown the smallest variance (V^-)min ~ (2j)^^'^ is obtainable at the time tmin ~ (2j)~^^'^. 
Based upon this, S0rensen et al. studied possible realization of the squeezing in ^'^Na 
atom BECs [Ij. More important, they proposed that the squeezing parameter can be 
used as a probe of many-body entanglement. In this paper, we investigate dynamical 
generation of the SSS in the generalized OAT model from arbitrary CSS. We find that 
the power rules change significantly even for ^o ~ 7^/2. 

At first, we determine the mean spin (J) = {{Jx), {Jy), {Jz)), where ( J^) = j cos(6'o), 
{Jx) = Re(J+), and (Jy) = lm(J+), with 

(J+) = j sin (6*0) exp [z(0o + 6t)] [cos (xt) + zcos {60) sin(xt)]^^~^ • (11) 
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It is convenient to rewrite Eq. (ITT]) as (7+) = rexp(i0), which yields (J^) = rcos0 and 
(Jy) = rsin0, as defined in Eq. ([1]). Therefore, we obtain 

r = J sin(^o) [1 - sin^ (^o) sin^ {xt)r'/\ (12) 
= 0o + 5t + (2j-l)(^(t), (13) 

where ip{t) = tan~^[cos(^o) tan(xt)] is dynamical phase. Note that in real calculations 
of the squeezing parameters, only cos(0) and sin(0) are needed and given by Eq. flA.2p 
and Eq. ( 1A.3[) . The explicit form of the phase or (/? is introduced to find out the 
roles of S and 0o in the squeezing. Obviously, r, if, and also i? = (r^ + {Jz)'^Y^'^ do not 
depend on them. 

To proceed, we calculate the expectation values {J+), and (J+(2J^ + 1)), which 
are relevant to the coefficients A, B, and C. The mean value (J|) reads 

( J|> = I sin^ (^o) + f cos' {9o) = \+3{3- 1/2) cos^ (^o) , (14) 

which, together with (J^) = jcos(^o); gives atom number variance (AA^^)^ = (AA^f,)^ = 
(AJ^)^ = (j'/2)sin^^o- Foi' ^/2, the variance (AJ^.)^ becomes narrow than that of 
the case 6^0 = 7r/2, which leads to relatively slow phase diffusion [l2l HS]. After some 
tedious calculations, we further obtain 

(4) = J (j - 1/2) sin2 (^o) exp [2z(0o + ^ 

X [cos {2xt) + i cos (^o) sin {2xt)f^~^ , (15) 

and 

( J+(2J, + 1)) = 2j (j - 1/2) sin (^o) exp [i{<P^ + 5t)] 
X [cos ixt) + i cos (^o) sin (x^)]^''"^ 
X [cos (6^0) cos ixt) + i sin (x^)] • (16) 

From Eq. (1A.4P -Eq. (lA.6p . one can find that the coefficients are fully determined by 
five quantities: sin^ (= r/R), cosO (= {Jz)/R), {Jl)e-'''^, and (J+(2J^ + l))e-^'^. 

We have shown that the first three terms are independent with 6 and 0o, which keeps 
true for the last two terms due to exp[i/c(0o + St)]e~^'"^ = exp[— zfc(2j — l)ip(t)] (with 
k = 1,2), where ip(t) does not depend on 6 and 0o- As a result, we get the conclusion 
that the detuning 6 and the azimuth angle 0o change the mean spin direction [see also 
Fig. [l](b) and (c)], but do not present any contribution to the squeezing. 

The squeezing parameters ^ and ( depend sensitively on the polar angle Oq of the 
initial CSS, as shown in Fig. [21 The most strongest squeezing can be obtained for 
^0 = ^/2, which corresponds to the initial CSS with equal atom population between 
the two internal states, i.e., (J^) = 0. From Eq. (fTTl) . we have (J+) = rexp{i(j)) 
with r = R = j cos^-'~^(x^) and (p = (po + St. From Eq. (fT^ -Eq. ( fT6l) . we further 
obtain (Jf) = j/2, {JDe'^'^ = 3(3 - 1/2) cos''^-\2xt), and (J+(2J, + l))e-'^ = 
13(23 — 1) cos^-'~^(xi) sin(xt). Substituting these results into Eq. (]A.4I) -Eq. (1A.6I) . 
we obtain the coefficients 

^=|(j-l/2)[l-cos^^--2(2xt)], (17) 



Spin squeezing in a generalized one-axis twisting model 



7 



B = 2j U - 1/2) cos^'~\xt) Mxt), (1^ 
and C = A + j . From Eq. ([5]) , we get the increased and the reduced variances 



I 
2 



1 + ^-^^ ± Va^ + 5^ 



(19) 



where the intermediate coefficients A = 1 — cos^-'~^(2xt) and B = 2 cos^-'^^(x^) sin(xt). 
One can find the variances are exactly the same with that of ideal OAT case [1], even 
for nonzero 6 and 0o. 




Figure 2. (Color online) Time evolution of the squeezing parameters (C^, for 
various of the initial CSS. From top to bottom: Oq = 7r/3 (squares, blue lines), 
0.98 xtt/2 (crosses, red lines), and tt/2 (empty circles, black lines). The arrows indicate 
the positions of the maximal-squeezing time imin for different ^o's. Other parameters: 
j = 30 (a), j = 2 X lO'' (b), and 5 = 0o = 0. 

Solid curves of Fig. [2] indicate the evolution of the normalized variance C,'^. The 
minimal value of the squeezing parameter, ^^^^ = 2j"^(Vl)min5 appears at the time 
tmin indicated by the arrows for different values of 6q. The smallest value of ^^j^ is 
obtained for the optimal initial state 6^0 = vr/2. For 7^ 7r/2 and large j (>> 1), the 
squeezing becomes worse than the optimal case. A closer look at the evolution of 
[= (j/ 1 (J) 1)^1^^] indicates that it is minimized before tmin [see empty circles of Fig. I^a)]. 
This is because different evolution rates of the variance V- and the mean spin (J). In 
addition, the minimal value Cmin is slightly larger than ^^^^ due to the decreased mean 
spin |(J)| < j. For large j case, however, the two squeezing parameters almost merge 
with each other in the short-time regime [see Fig. [2]|^b)]. As a result, one can assume 
that ^j^jj^ obeys the same power rule with C,^^^ [M], and is determined by that of (y_)mm- 

4. Power rules of the strongest squeezing and its time scale 

As shown by the red lines of Fig. [2]^b), both ^^j^ and tmin change significantly in 
comparison with the idea case (i.e., 6*0 = 7r/2). As a result, it is necessary to determine 
power rules of the variance (\^_)min and the time tmin for Oq vr/2. In this section. 
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we calculate analytically the power rules by using standard treatments of Ref. pQ. We 
will focus on a small departure of from 7r/2 due to the fact that a relatively small 
population imbalance between two internal states favors the one-axis twisting effect. 

4-1. Ideal OAT case with 6q = ir /2 

In the short-time limit {xt << 1) and large particle number {j >> 1), the increased 
and reduced variances Eq. ( ITQll can be approximated as [1]: 

V,^l(,o$),V.^l{^^l^, (20) 



where ao = jxt > 1 and /?o = iix't)'^ << 1- Eq. (j20l) is the key point to obtain the 
strongest squeezing and its time scale t^^-^. Previously, the time tmin was obtained 
by comparing the second term of with that of the first one [1]. Here, we solve tmin 
via minimizing V- with respect to t, i.e.. 



dt^ 



0, (21) 



which yields power rule of the maximal-squeezing time: 

Xtu.,n^^''\2j)-^l\ (22) 

Inserting x^min into Eq. fl20l) . we further obtain the reduced variance as 

2,xV3 



[V^Un ^ 8 (^fj ' (23) 

and also, the smallest squeezing parameter ^^j^^ = 2j^^(V_)inin ~ \(^)~'^^^- Power 
exponents of Eq. fl22l) and Eq. fl23l) are consistent with Ref. [1], but different in the 
coefficients. As shown by the black solid lines of Fig. [3l the revised results fit very well 
with their numerical results (empty circles). 

4-2. Small departure case with 6q ~ it/2 

The power rules, Eq. ( 1221) and Eq. (l23l) . are valid only for 6q = 7r/2. Now, we generalize 
them for 6q ^ 7r/2 case. To obtain the approximated expressions of the variances as Eq. 
( l20l) . we calculate short-time solutions of (^+), and {J+{2Jz + 1)). 

In the short-time limit (x^ << 1), the dynamical phase ^p{t) = 
tan~^[cos(^o) tan(xt)] 2^ xtcos(6'o), and Eq. (ITT]) can be approximated as 

(J+)^jsin(^o)e*V'^, (24) 

where [3 = (3q sin^(6'o) = jix^Y sin^(6'o), and </> ~ (f)o + St+2jxt cos(6'o). We have assumed 
that particle number is large enough so 2j — 1 ~ 2j. The length of the correlation 
reads r = — jsin(^o)e~^, which indicates that phase coherence of the two-mode 

BEC decays exponentially (i.e., phase diffusion [11]) with the coherence time scaled as 



Spin squeezing in a generalized one-axis twisting model 



9 



Xtd = sin ^(6'o)j ^^"^ [121113]. Similarly, short-time solutions of Eq. (IT^ and Eq. (HM 
can be written approximately as 

{jl)c^j{j-l/2)sm^9o)e''h-'^, (25) 

and 

(J+(2J, + 1)) ~ j(2j - 1) sin^olcos^o + ^Xt)e"^e-^, (26) 

where the factor cos^o can not be neglected since it is comparable with xt- In fact, 
it is the presence of cos 6^0 that leads to significant change of tmin and (\^_)min even for 
60 ~ vr/2. 




Figure 3. (Color online) The maximal-squeezing time tniin (a), and the smallest 
variance {V-)min (b) as a function of j (= N/2) for Oq = n/3 (squares, blue hues), 
0.98 X 7r/2 (crosses, red lines), and tt/2 (empty circles, black lines). Solid lines are 
predicted by Eq. (|28p and Eq. (P^)) . Dashed lines are plotted to guide the eyes, and 
are given by i(2j)~^/^ (a) and {2j)^^^ (b), respectively. Other parameters: S ~ (I)q — 0, 
and iniin is in unit of x~^- 



To simplify the calculations, we make further approximations to the angles of Eq. 
(lA.4p -Eq. (1A.6P : sin^ = r/R ~ sin6'o and cos6' = {Jz)/R — cos^qj where Oq is 
polar angle of the initial CSS. This approximation is equivalent with r ~ j sin(6'o), i.e., 
neglecting the the diffusion within the squeezing time due to td > tmin- Now, we expand 
the coefficients A, B, and C in terms of /3. In calculating the increased variance, we only 
keep the lowest order of P, and get V+ ^ |(4q;^), where a = aosin^^o = JX^sin^^o- 
Next, we solve power series of AV^V^ up to the third order of /?, from which we obtain 
the reduced variance as 

2/52 



2 



+ ^ (l + Qjsin^^o cos'^o) 



(27) 



_4a2 ■ 3 

where the j-dependent additional term gives significant contribution to the squeezing 
for 6*0 7^ 7r/2. By minimizing with respect to t, we obtain power rule of the time as 

_ 3^/6(2jsin2 60)-^/-' 

(l + 9jsin2 00 cos2 ^o) 
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and that of the decreased variance: 



2j 



3 sin 9q 



(1 + 9jsin2 ^0 cos^ ^o) 



1/3 



(29) 



For 6*0 = 7r/2, our results reduce to the ideal OAT case, i.e., Eq. 0221) and Eq. 
while for 6*0 7^ 7r/2 and large j, Eq. (!28l) and Eq. (l29l) predict that the power rules 
change to 

Xtmin ~ i2j)-'/', {V.U, ~ i2jf\ (30) 

which are confirmed by numerical simulations. To see this more clearly, let us focus on 
red lines of Fig. [3l For 60 ~ it/ 2 and small j, both the time x^min and the variance 
(^-)min follow the Same rule with the ^0 = case. With the increase of j, however, 
the red line (the crosses) of Fig. decreases faster than the ideal OAT case [see also 
Fig. [2](b)]. The change of the power rule is shown more clearly in Fig. [3](b). 

In Fig. m we show the dependence of tmin and (l^_)min on 9q for a fixed value j. 
It was show that both tmin and (l^-)mm are symmetrical with respect to 9o = 7i/2. 
The most strongest squeezing [i.e., the smallest value of ^^j^] occurs for the optimal 
initial state 6^0 = 7r/2. Our analytic results, Eq. (j28l) and Eq. ( j29ll . agree quite 
well with numerical simulations except 6*0 = or tt. In this case, the state vector 
|\l/(t)) = exp[— z(xj^ ± Sj)t]\j, ±j), which is the CSS with the variances (V^+) = (V^) = 
j/2 and = 1. However, Eq. ( l29i) diverges as 6'o — or vr, inconsistent with 
the real situation. Eq. (128!) gives relatively good estimate of the maximal-squeezing 
time. As shown in Fig. Hl^a), tmin decreases monotonically in the small departure regime 
1^0 ~ 7r/2| < 0.277r/2, which implies that the maximal squeezing occurs more and more 
earher [see also Fig. [2](b)]. Out of the regime, tmin increases with the departure of ^0, 
and goes infinity as 6'o — or tt. 



(a) 





00 



Figure 4. The time tmin (a) and the normahzed variance ^"^^^ = 2(y_)niin/j (b) as 
a function of (in unit of 7r/2) for j = 10^. Empty circles are given by numerical 
simulations, and solid lines are predicted by Eq. ([28)) and Eq. ((29l) . respectively. 
Vertical grid lines in (a) denote 9q = 0.73 x 7r/2 and 6*0 = 1-27 x 7r/2, separating 
monotonic regimes of tmin- Other parameters are the same with Fig. 3. 
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4-3. Dissipation effect due to atomic decay 

So far, we have neglected the effects of dissipation on the spin squeezing, such as particle 
losses and center-of-motion of the atoms in the BECs [H HOI SI]. For the squeezing 
generated in atomic ensemble, the dominant dissipation source is atomic decay due to 
spontaneous emission [32], which can be described by the master equation [26] : 

^ = z[p, if] + I (2 J_pJ+ - J+ J„p - pj+ J_) , (31) 

where p is the density operator, and 7 is the decay rate of the atoms. In the basis of 
\j, m), the elements pm,n = {j, '^IpIj, n) could be solved numerically by using the Runge- 
Kutta routine [26] . In real calculations of the squeezing parameters, only 6j elements 
hke pm,m, Pm,m+i, and Pm,m+2 are needed. 

In Fig. [5l we plot time evolution of for small decay rate, e.g., 7/x = 0.01 and 
0.1. Such a small dissipation can be realized by increasing x, which in turn leads to the 
preparation of the SSS within the lifetime of the atoms 7"^ [32]. For relatively small 
decay rate 7/x = 0.01 (red curves), both the maximal squeezing and its time scale 
change slightly in comparison with 7 = case. The initial state with = 7r/2 looks 
more sensitive to atomic decay than 9q = 0.8 x n/2 case. From the blue dotted lines of 
Fig. [5](b), we find that even for 7/x = 0.1, a considerable squeezing with (and also 
C^)~ 0.22 could be reached in an ensemble of 200 atoms, which occurs at a time scale 
given by Eq. 



(a) (b) 




Figure 5. (Color online) Time evolution of for 0^ = 7r/2 (a) and 6*0 = 0.8 x 7r/2 
(b). From top to bottom, the decay rates 7 = O.lx (blue dotted line), 7 = O.Olx (red 
solid line), and 7 = (black dashed line). Vertical and horizontal grid lines denote, 
respectively, ^min and f^j^ for 7 = case. Other parameters: j = 100, (5 = 0o = 0. 



5. Conclusion 

In summary, we have presented general formulas to study spin squeezing in spin- 1/2 
system. Instead of six fluctuation parameters as Refs. [IQ], HI], only five parameters, 
i.e., (J2), {Jl), (</+), and {J+{2Jz + 1)) are needed to determine the mean spin 

and the squeezing parameters. 
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Spin squeezing of a generalized one-axis twisting model is investigated for arbitrary 
coherent spin state 16*0, (po). We show explicitly that 6*0 = vr/2 is the optimal initial state 
to obtain the minimum value of the variance {V-)mm — |(2j/3)^/^, which takes place at 
the time scaled as x^min — 3^/^(2j)~^/^. The detuning S and the azimuth angle 0o alert 
the mean spin's direction, but give vanishing contribution to the squeezing parameters. 
As the main result of our paper, we calculate analytically the dependence of the variance 
(V_)min and the time tmin on the polar angle 6*0, as Eq. fl28l) and Eq. fl29|) . respectively. 
What may be a little surprising is that even for a small departure of Oq from 7r/2, 
the power rules become (V^_)min ~ i'^jY^^ and x^min ~ (2j)~^^^, deviating from the 
ideal case. The power rule, Eq. ( l28l) . is robust against atomic decay for 7 < O.lx- 
Our results show that spin squeezed state in the OAT model depends sensitively on the 
initial state and the interaction time. A straightforward way to overcome these stringent 
requirements is still an open question. 
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Appendix A. The coefficients A, B, and C 

In Eq. ([T]), we have defined three orthogonal unit vectors n^ {i = 1,2,3), which are 
valid for any spin state l^'). The angles 6 and are determined by the mean spin 

sin9 = l cos9=i^, (A.l) 

c„s^=<|>=M^. (A.2) 

_ (i> ^ (A.3, 

where the length of the mean spin R = |(J)| = {{Jx)"^ + (Jy)"^ + {Jz)"^)^^"^ and 
r = |(J+)| = ((4.)2 + {Jy)Y^^ = RsinO. From Eqs. flATTjl -fm. it is easy to verify 
that the unit vector = (sin ^ cos 0, sin 6* sin 0, cos ^) = i?~^(J), i.e., the mean spin 
(J) is parallel with the unit vector n^. Moreover, one can prove the expectation value 
(-^ni) = -(Jx)sin0+ {Jy)cos(p = -{Jx){Jy)/r + {Jy){Jx)/r = 0, (Jna) = 0, and 

By using the above results, one can solve explicit expressions of the coefficients: 
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A = {Jl^ - JlJ, B = (Jni^2 + Jn^Jni), and C = {Jl^ + J2^), yielding 

2A = sin^ e [jU + 1) - 3 ( J2)] - (1 + cos^ 6) Re [( ) e'^^^j 

+ sin(2e) Re [( J+(2J, + 1)) e"*"^] , (A.4) 

S = - cos(e) Im [{Jl) e-^'^] + sm{e) Im [( J+(2J, + 1)) e"*^] , (A.5) 

C + A = jij + 1)- {J',) - Re [(4) , (A.6) 

where we have used the relations: (Jj + Jy) = j{j + 1) — (Jf), {J^ — Jy) = Re(J^), 
{J^Jy + JyJ^) = Im(J^), (J^J^ + J^J^) = Re(J+(2J^ + 1)), and {JyJ^ + J^J^) = 
Im(J+(2J^ + 1)). Substituting the coefficients into Eq. ([5]), we obtain the variances 
V± and the squeezing parameters and 

Now let us calculate the coefficients for any CSS \0,(j)). The mean values 
{Jz) = jcos(6') and {J^) = j sin (6') e**^, which yields (J^) = jsm{9)cos(j) and 
(^j/) = J sin (6*) sin</). For the CSS, the length of the mean spin R = j. These results 
can be directly obtained from Eq. (fTTi) by taking time t = 0. Similarly, from Eqs. 
([I4D-([T6D, we further obtain {Jl)e-^''^ = j{j - 1/2) sin2(^) and (J+(2J^ + l))e-'^ = 
j{2j — 1) sin(^) cos(6'), from which we obtain immediately the coefficient B = 0. 
Substituting these results into Eq. flA.4l) -Eq. flA.6p . we also get ^ = 0, and C = j. 
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